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Abstract 



^SJ ' The CP-violating asymmetries in the exclusive decay B — > K*rl^ (/ = 

^ . e^iJL,T) are predicted to be exceedingly small in the standard model (SM), 

thereby offering an opportunity to assess various new-physics scenarios. We 
CN ' derive quantitative predictions for various integrated observables m. B ^ 

^ I K* ijl'^ (JL~ decay in the presence of physics beyond the SM with additional CP 

O ' phases and an extended operator basis. In particular, a model-independent 



analysis of CP asymmetries that require the presence of unitarity phases, in 
addition to CP violation, is performed. We find that in the low dimuon in- 
variant mass region 2m ^ ^ M^+^- < Mj^^, the CP asymmetries are highly 
suppressed by small dynamical phases, assuming that new physics is unlikely 
' to significantly alter the Wilson coefficients of the operators governing two- 

body hadronic B decays. Taking into account current experimental data on 
^ , the measured 6 — > 57 rate and the upper limit on B{B^ K*^fi^fj,^), CP- 

^ I violating effects of a few per cent are estimated, even in the presence of new 

physics with CP phases of 0{1). By contrast, in the high dimuon invariant 
mass region M^/ < M^+^~ ^ {Mb — Mk*) significant CP-violating effects 
are possible. Given a branching ratio of 1.8 x 10^^, the CP asymmetries can 
be quite substantial (~ 20% or more), and thus may serve as a means of 
discovering physics transcending the SM. 
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I. INTRODUCTION 



CP violation has been observed so far only in the neutral kaon system. Within the stan- 
dard model (SM), the experimental results on indirect (ex) and direct (e'/ex) CP violation 
can be explained by the complex phase of the Cabibbo-Kobayashi-Maskawa (CKM) matrix 
if one takes into account the large theoretical uncertainties associated with the hadronic 
matrix elements that enter the analysis of e'/ex [Q- It therefore remains an open question 
whether the CKM mechanism of CP violation can account for the new experimental result 
on e'/ex- 

A great deal of effort is given to the study of CP violation in the B system, which will 
provide invaluable information on the pattern of CP violation and open up the possibility 
to look for new physics. In this paper we are concerned with the exclusive decay B — *■ 
K*l^l~ which is of special interest because (i) it probes the underlying effective Hamiltonian 
describing flavour-changing neutral current (FCNC) processes in B decay; (ii) the analysis 
of CP-violating effects in decays governed by 6 ^ sl^l~ may offer a deeper insight into 
the mechanism of CP violation since the SM prediction for CP asymmetries is extremely 
small, typically < 10^^ 0; and (iii) the process B — > K*fi^fi~ is a very promising decay 
mode since it is likely to be observed in the next round of B physics experiments. The most 
stringent limit has been set by CDF of B{B^ K*°i2+fi-) < 4.0 x lO'^ at the 90% C.L. ^ 
to be compared with the SM prediction of (1.9 ± 0.7) x lO"*^ [|]. 

The object of the present work is to explore the possibility of sizable CP asymmetries 
in i? ^ K*l^l~ decay, whose observation would clearly indicate the presence of physics 
beyond the SM. We perform a largely model-independent analysis by considering a new- 
physics scenario with additional CP phases and an extended operator basis, taking into 
account existing experimental data on B ^ and the upper bound on B^ K*'^fi^fi~ . 

Our paper is organized as follows. Section ^ contains the SM operator basis and short- 
distance matrix element for the process B K*l^l~ . A formalism for dealing with real cc 
intermediate states such as J/il),ip' entering via the decay chain B K*Vcc — ^ K*l~^l~, is 



described. In Sec. |T|, we briefly discuss possible extensions of the SM including those with 
an extended operator basis. Section jTV] is concerned with the parametrization of the hadronic 
matrix elements and gives the differential decay spectrum as well as the forward-backward 
asymmetry. The corresponding CP- violating observables are discussed in Sec. |V]. Numerical 
estimates for integrated observables in 5 — »• K*n'^n~ decay in the presence of new physics 



are given in Sec. VI . Particular attention is paid to the CP- violating partial-rate asymmetry 



and the CP-violating effect related to the angular distribution of in i? and B decays. 



Our conclusions are contained in Sec. VII. 



II. EFFECTIVE HAMILTONIAN 

A. Short-distance contributions 

The effective Hamiltonian for the decay B K*l^l~ in the standard model (SM) is 
5iven by 
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(2.1) 

where we have used the unitarity of the CKM matrix, = Vub^usl'^tb^tsi ^^^1 the operator 
basis is defined as follows: 

01 = {Sal^.PLCp){cpl^PLha): 

02 = is^l^PLCa)ic^l^PLbl3), 

= {s^^^PLUa){upYPLhp), 

q=u,d,s,c,b 
q=u,d,s,c,b 
q=u,d,s,c,b 

0, = {-s^l,PLhp) Yl (K^^'Pr^c.), 



q=u,d,s,c,b 



O7 



167r2 

9s 
167r2 

167r2 
e2 



Sa(^,,u{mbPR + msPL)haF^\ 

Sal^PhhahtJ-, 
-,Sa1^PLhJ'^^,'^d, 



(2.2) 



where a, (3 are colour indices, a labels the SU(3) generators, and Pl,b. = (lT75)/2. Evolution 
of the Wilson coefficients Cj(/i) in Eq. ( p.l|) from the weak scale /x = Mw down to the low 
energy scale fi = rrih by means of the renormalization group equations (RGE's) then leads 
to the QCD-corrected matrix element in next-to- leading logarithmic approximation 



M = 



tb^ts 



cf-cio){K*{k)\s^,PLb\B{p)) 



2cf 



{K*{k)\-si(j^,q' {rribPR + ttIsPl) h\B{p)) 



I^PlI + (cio ^ -cio)irPRl }■ (2.3) 



Here q = p — k,s = q'^ is the invariant mass of the lepton pair, and the effective Wilson 
coefficient cf^ has the form 



cf = c, + Y{s), 



(2.4) 



with 
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TABLE I. Numerical values of the Wilson coefficients ci, . . . , cio at /i = mf, within the SM. 



Cl 


C2 


C3 


C4 


C5 


C6 


„eff 




ClO 


-0.249 


+1.108 


+0.011 


-0.026 


+0.007 


-0.031 


-0.314 


+4.216 


-4.582 



= g{mc, s)(3ci + C2 + Scs + C4 + Scs + ce) + \u[gimc, s) - g{mu, s)](3ci + C2) 

1 1 

- -g{m„ s)(c3 + 3C4) - -g{mb, s)(4c3 + 4c4 + 805 + Ce) 

+ ^(3C3 + C4 + 3C5 + C6) + ---, (2.5) 

where the ellipsis represents the order «<, correction to the matrix element of the operator Og, 
which can be regarded as a contribution to the form factors [0 , and hence will be omitted in 
the calculations that follow. Table |I| summarizes our results for the Wilson coefficients Cj(mf,). 
Observe that cf^, Cg, and Cio are real in the framework of the SM. The function g{mi, s) in 
the above formula arises from the one- loop contributions of the four-quark operators Oi-Oq, 
and is given by (at /i = rrib) 
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g{mi,s) = --\n{mi/mb) + ^ + gV^ ~ -(2 + V|l - 



e(i-y. 



ln(i±^l-.. 



QiVi - 1)2 arctan } , (2.6) 



with Hi = Amfls. This expression reduces to 



g{^.s) = ^-Un{s/ml) + '^i7,, (2.7) 

in the limit mj ^ 0. A few remarks are in order here. 

(i) The Wilson coefficient Cg^, Eq. ( |2.4| ), has absorptive parts for s > 4m\ and s > 4:ml, 
and thus contains dynamical (unitarity) phases. As will become clear, these are prerequisites, 
besides a CP-violating phase, for observing CP asymmetries in partial rates. Since the 
remaining coefficients cf^ and Cio do not contain any strong phases, the unitarity phases 
below the cc threshold are generated by light quark contributions, whereas for s > Arn^ they 
arise mainly from cc intermediate states. 

(ii) Within the SM, CP violation in b sl~^l~ transition is caused by the ratio of CKM 
factors appearing in cf^, namely A„ ~ (A = Vus — 0.22), which is further reduced by a 
factor of order (3ci + C2)/cg ~ 0.085. As a result, the effective Hamiltonian for h —>■ sl^l~ 
essentially involves only one independent CKM factor VJ^VJ^, so that CP violation in this 
channel is unobservably small. Thus, the numerical effect of A^ in Eq. ( p.5| ) is negligible for 
decays based on the transition h sl^l~ . 



B. Resonance contributions 

In addition to the short-distance contributions discussed so far, there are possible quark 
antiquark resonant intermediate states like (p, J/tp,ilj' etc. Since the s-quark contributions 



4 



in Eq. ( p.5|) are suppressed by small Wilson coefficients, and terms proportional to may 
be dropped in the case of 6 — s transition, we are left with the J /ip family. 

Following the procedure in Ref. p[, we implement the charmonium resonances utilizing 
e^e~ annihilation data. The absorptive part of the one-loop function g, Eq. ( p.6|) , can be 
related to the experimentally accessible quantity 

R[s) = crtot(e^e^ — > hadrons)/cr(e^e^ — > yU^/i^) (2.8) 

by virtue of the optical theorem. Specifically, it is found that 

lmg{m,,s) = ^R''/^{s), (2.9) 

with R'^/^{s) = i?^ont(s) + Rri^{s), whereas the dispersive part Reg{mc, s) may be obtained 
via a once-subtracted dispersion relation 



s 



oo 



g{m,, s) = g{m,, 0) + J / J^'^'^^'K ds', e +0, (2.10) 

^ JiM'^ ^ \^ — s — le) 

with (yf(mc,0) = — 8/9 ln(mc/mfc) — 4/9. The continuum contributions, -Rcont; can be deter- 
mined using experimental data from Ref. 0, while the narrow resonances are well described 
by a relativistic Breit-Wigner distribution P,pil[ 



with the properties of the vector mesons summarized in Ref. |TT[ 



To account for experimental data on direct J/ip production via the relation 

B{B K*Vcc K*l+l-) = B{B K*Vcc)B{Vc-c /+r), (2.12) 

where Vcs = J/'ip, "ip', ■ ■ ■, one usually modifies the Breit-Wigner distribution in Eq. (|2.11|) by 
introducing an ad hoc factor Ky 0, 12|. This suggests that the factorization ansatz inherent 



in the approaches that have been advocated to incorporate resonance effects in 6 ^ sl^l 
(see, e.g., Ref. |]13[) is inadequate for two-body hadronic decays of B mesons. Using the set 



of form factors that will be discussed below, one finds kj/^ = 1.7, k^' = 2.4, whereas for the 
remaining resonances we shall take (as in Ref. 0]) k = 2. In Fig. |1], we show the real and 
imaginary parts of g{mc, s) based on Eq. (|2.10|) , as a function of s = s/M"^. 



III. EFFECTIVE HAMILTONIAN AND NEW PHYSICS 

The new-physics contributions to the decay mode h —>■ sl^l~ can manifest itself in two 
distinct ways: (i) The absolute values and phases of the Wilson coefficients at the electroweak 
scale are modified, (ii) New operators in addition to the ones defined in Eq. ( p.2| ) arise. We 
consider them in turn. 
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FIG. 1. The predicted s dependence of g{mc,s) using experimental data on e~^e~ — > hadrons 
via a dispersion relation, Eq. ( 2.10| ), including the effects of cc resonances (solid curve). Note that 
Im g{mc, s) 7^ is entirely due to unitarity phases. Also shown is the result of the perturbative 
calculation according to Eq. (^) (dashed curve). 



A. Wilson coefficients and new physics 



The non-standard contributions can be divided into three main groups |]I4[: First, models 
with tree-level contributions to the four-quark operators such as supersymmetry (SUSY) 
with broken R-parity or models with Z-mediated flavour-changing neutral currents. Second, 
models with contributions at one-loop level with new particles running in the loop like 
SUSY with conserved R-parity or models with four quark generations. Third, models with 
no significant effect on the above Wilson coefficients including multi-Higgs-doublet models 
with natural flavour conservation and left-right symmetric models. 

Let us begin with the Wilson coefficients C3, . . . , of the QCD penguin operators. It fol- 
lows from the RGE analysis that their numerical values at /i = rrib are essentially determined 
by the Wilson coefficient C2 of the four-quark operator O2 at the electroweak scale. Thus, 
in order to have considerable deviations from the SM predictions for the coefficients cs-cg, 
one needs large new-physics contributions to the short- distance coefficient C2, which in turn 
would affect the theoretical branching ratio for two-body non-leptonic B decays. Recent 



studies |T5| suggest, however, that the short- distance coefficients of the SM can account for 
existing data if one allows departures from the naive factorization prescription. 

For the numerical analysis here, we adopt the SM values for the Wilson coefficients 
Ci, . . . , Cg summarized in Table |. Furthermore, it is convenient for later use to parametrize 
the new-physics contributions to the remaining coefficients cf^, c^, and Cio by the following 
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ratios (defined at the scale /i = rrib): 

R, = cf/cf'^^ = \R7\e'^\ (3.1a) 

Rg=Cg/4^ =\Rg\e'^\ (3.1b) 

Rio = = \Rio\e"^''. (3.1c) 

where we recall Eq. (|2.4|) , and 

c, = cS^^ + cf^-, (3.2) 

with cf^ shown in Table |. Note that the Wilson coefficient Cio in the above does not depend 
on the renormalization scale, and thus Cio = ciQ^Myy). 

B. Extended operator basis 

It is conceivable that physics beyond the SM induces new operator structures containing 
scalar, pseudoscalar, and tensor interactions, in addition to the SM operator basis, Eq. (^.21). 
The Wilson coefficients of the scalar operators 

Of = SaPLbJl, 

of = s^PnbJl, 

Ofo = So^PRhaHl, (3.3) 

involve the lepton mass in most extensions of the SM, and hence will give only small contri- 
butions in the case of / = e or /i. Moreover, possible tensor-type operators sa^^hla^^l and 
sia^yhlaapl^^^""^ can be safely neglected since their numerical contributions have been found 
to be small |]TB[. 

This leaves an extended operator basis which consists of the SM operators and the 
opposite-chirality operators p!6|,p!7[ 



^7 = -r^s^a^,{m,PL + m,P«)6,F'^^ 

lOTT"^ 

O; = -^7s^T^f,a^,{m^PL + m,Pfi)6;3G"^'^^ 

Og = y^s^i'^PrKIi^i. 

^io = T7r^s^l^PRbc.ll^.l^l. (3.4) 
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TABLE II. LCSR predictions for the B ^ K* form factors with f{s) = /(O)exp(cis + C2s'^), 
s = s/M-^ Q. Recall that is given in terms of Ai and A2 via Eq. (4.3). 





V 


Ai 


A2 


Ao 


Ti 


T2 


n 


/(O) 


0.457 


0.337 


0.282 


0.471 


0.379 


0.379 


0.260 


Cl 


1.482 


0.602 


1.172 


1.505 


1.519 


0.517 


1.129 


C2 


1.015 


0.258 


0.567 


0.710 


1.030 


0.426 


1.128 



IV. DECAY DISTRIBUTION FOR B K*l+l- 
A. Form factors 



The hadronic matrix elements appearing in Eq. ( |2.3[ ) can be expressed in terms of s- 
dependent form factors (recall s = g^); namely 

{K*{k)\s^,PL,Rb\B{p)) = ^e,.«^e>"g^^^^-^^t^ T ^{e;(M^ + MK^iis) 

- (6- ■ q){2p - - ^{e* ■ q^sis) - Ao(.)]g,| , (4.1) 

with the convention €0123 = Q = P — k, and 



{K*{k)\7sia,,q''Pn,Lb\B{p)) = -ze^,„^e>"g/'Ti(s) ± 1<] [e;(M| - 



2 



-(e*-g)(2p-g)jT2(s) + (e*-g) 



Tsis) , (4.2) 



where Ti(0) = T2(0), and is the K* polarization vector. The form factor A3 can be 
written in terms of Ai and A2, i.e. 



with the relation ^3(0) = Ao(0). The terms proportional to in Eqs. (f4.1|) and ( [4.2|) 
do not contribute to the differential decay rate when the final leptons are massless. For 
the results presented below, we adopt the B —>■ K* form factors of Ref. which have 
been obtained using hght cone sum rule (LCSR) results, and are displayed in Table f|. 
Throughout our discussion, we assume the above form factors to be real, in the absence of 
final-state interactions. 



B. Differential decay spectrum 

Squaring the matrix element (^.31), summing over spins, and introducing the shorthand 



notation 

,2 , 1,2 , „2 



A(a, 5, c) = a^ + ¥ + c- 2{ab + bc + ac), (4.4) 
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Mi = M,/Mb, rhi = nii/MB, s = s/M\ 



(4.5) 



(4.6) 



the spectrum of S ^ K*l^l decay with respect to s and 9^ the angle between / and the 



outgoing hadron in the dilepton centre-of-mass system, is |T9| 



The quantities A, B, C are defined as follows: 
2X2 r 



[A{s) + B{s) cos 01 + C{s) cos^^/]. (4.7) 



Ais) 



sMhHs) + - ( 1 + ) Us) + X^m + U{s) 



+ 2mf/(S), (4.8) 



B{s) = 8X 



4m? 
1 ^ 



{'^10 [''9 ^^xAy Cj [AxBy + AyBx)]} 1 



(4.9) 



Cis) 



2X2 

Mi, 



4m2 



sMl,m - -/2(5) - XV3(S) - U{S) 



with the auxiliary functions 

I{s)=AX'Ms) + Ms) + Ms), 



(4.10) 



(4.11) 



/i(5) = (|cfP+|cior)A2 + 



4\cf\[^, 4Re{cfcf*) 



AxBxi 



(4.12) 



/2(5) = W). Ms) = Ms)^_ 



>2' 



(4.13) 



his) = - s - M'^.) 



(|cf r + + ^^ByB^ - [cfcf*) ^^^^^ ^ ^^^^^ 



(4.14) 



Ms) = IcioP 



^X^A^ - 3Al + ^{[2(1 + Mi.) - + 2Ay}A, 
4X 



+ -^{M^.(A3 - Ao) - [(1 - MiO^ + ^,]}(^ - A 



sM, 



K* 



(4.15) 



In these equations, 
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Vis) 



1 + Mk. 



Ay = {1 + Mk')Ai{s), a. 



A,(s) 



1 + 



(4.16) 



= -Ti{s){mb + rhs), By = -(1 - M^,)T2{s){mb - rhs 



(4.17) 



B. 



Us) + 



1-M|. 



{rhb - rhs), 



(4.18) 



with the form factors hsted in Table y. 

A further observable of interest is the forward-backward (FB) asymmetry of /~, defined 

as 



Afb{s) 



acos^i— - 

as a cos Ui 



° dV 
d cos 61—— - 

_i as a cos ai 



dcosdi-— - 

as a COS Ui j 



° dV 
d cos 61 - 



(4.19) 



' dsdcosdi 



and we obtain 

^fb(s) 



,0^, A 4mf Re{c-o[cfM.A,-cf(A,.i?, + A,i?,.)]} 
s [?>A{s)+C{s)] 



(4.20) 



Note that the Wilson coefficients in the above expressions are defined through Eq. (|3.2| ). The 
decay distributions in the presence of new physics with additional operators are discussed 



further in Sec. 0. We now proceed to a discussion of CP-violating observables in the process 
B K*l+l-. 



V. CP- VIOLATING OBSERVABLES 

Suppose the decay amplitude for B F has the general form 

A{B ^F) = e'^^Aie'^^ + e'^''A2e'^\ (5.1) 

where Ai^2 are real matrix elements, 5i and 0j are the strong phases (CP-conserving) and 
weak phases (CP-violating) respectively. Using CPT invariance, which requires that the 
total decay rate for particle and antiparticle be equal, the decay amplitude for the conjugate 
process takes the form 

'A{B ^ F) = e'^'^^Aie'^^ + e-'*^A2e'^\ (5.2) 
giving rise to the CP asymmetry 



I^P - I^P _ -2rsin0sin(5 
|^|2 + \A\^ ~ 1 + 2r cos0cos5 + r2 



with r = A2/A1, (f) = 01 — 02; and 5 = 5i — 82- Notice that in the limit r <^ 1 the asymmetry 
is approximately Acp ~ — 2r sin sin 5. Inspection of Eq. (|5.3|) reveals that a non-zero 
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partial-rate asymmetry requires CP violation (0 7^ 0) and the presence of dynamical phases 
{S 7^ 0), the latter being provided by the one-loop function g{mi,s) present in the Wilson 
coefficient cf [Eq. (^]. 

To determine the impact of the strong phases on the CP asymmetries, it is useful to 
separate those contributions to the decay amplitude that generate absorptive parts. To this 
end, we choose 5i = 02 = in Eq. (|5.1|) and require that A^e^^'^ vanishes when Yi^s) 0. 
Moreover, as discussed in Sec. fTT|, the corrections to the short-distance coefficients of the SM 
multiplying the absorptive parts of the decay amplitude are not expected to be large in many 
extensions of the SM, nor required by current data on two-body hadronic B decays ||15| . 
Consequently, the strong phases in 6 ^ sl^l~ decay are essentially unaffected by possible 
new interactions transcending the SM, so that the non-standard effects on CP asymmetries 
can be described by the two parameters r and 0. 

In order to get a quantitative idea of the magnitude of unitarity phases in 6 — >■ sl'^l~ 
transition, we plot in Fig. |^ the parameter Ry = — Im [y(s)]/|y(s)| for two different regions 
of the dilepton invariant mass, namely 

4m^^s<: {Mj/^ - e^ut)', (M^' + e^t)' ^ ^ ^ (1 - Mk*)\ (5.4) 

which we refer to as the low-s and high-s region respectively.^] As we alluded to earlier, 
in the low-s region the strong phase is suppressed by small Wilson coefficients of the QCD 
penguin operators, whereas in the high-s region it can be large. 

Figure |^ shows the dependence of A^-p on r for different choices of the weak phase sin 
in the low-s and high-s region. Observe that the asymmetry is maximized when r = 1 (i.e. 
the two interfering amplitudes are of comparable size) but is strongly suppressed if either 
r ^ 1 or r ^ 1. 

Using the two-dimensional decay distribution dT/ds d cos 61 derived in the preceding sec- 
tion, we may define the following average CP asymmetries: 

ds / a cosy;—— — / ds d cos Ur 



/ aD \ _ Jso Jd dsdcosOi / \ _ Jsp Js dsdcosOi 



«i dT sum ^ l-si /• + ! dVsnm 

ds / dcos6i /"""^ / ds dcosOi /""^ 

so J-i dsd cos 9i Jg^ 7_i dsd cos 9^ 



where 



L 



-1 pO 

^ I T , (5.6) 



D,S Jo J-1 



Fsum = r{B ^ K*l+l-) + r{B ^ K*l+l-), (5.7) 
Fdiff = T{B K*l+l-) - T{B K*l+l-). (5.8) 



'^We use ecut = 0-2 GeV/M^ and e^ut = 0.1 GeV/M^ in the case of in the final state. 
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0.06 



(a) sin S = —0.1 



sin (j> = 1 
sin (/> = 0.6 
sm6 = 0.2 




0.3 



(b) sin 5= -0.5 



sin 0=1 
sin (j> = 0.6 
sin 6 = 0.2 



4 6 

r 



FIG. 3. CP- violating asymmetry ^CP) Eq. ( ^.3[ ), as a function of r for (a) the low-s region with 
sin (5 = —0.1 and (b) the high-s region with sin (5 = —0.5. 
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Notice that the asymmetry A^p is a CP- violating effect in the angular distribution of /~ 
(or equivalently the Dalitz-plot distribution) in B and B decays while A^p represents the 
asymmetry in the partial widths of these decays. The special feature of the former is that 
it can be determined even for an untagged equal mixture of B and B events, i.e. without 
flavour identification. We emphasize that both asymmetries are odd under CP but even 
under 'naive' T, and thus vanish in the limit that there are no strong phases. [| 

VI. NUMERICAL ANALYSIS 

In the remainder of this paper, we wish to focus on the B — > K*n~^fi^ mode. We start 
this section with a brief discussion of the experimental constraints that we shall be using in 
our subsequent calculations. 

A. Experimental constraints 

Constraints on the parameters Ri, Eqs. ( ^.11 ), are provided by the following experimental 
results: (i) The CDF collaboration has recently obtained the upper bound 

13{B^ K*^fi+fi-) < 4.0 X 10"^ (6.1) 

at 90% C.L. 01 . (ii) The measurement of the inclusive branching ratio B{B X^^y) yields 

m 

2.0 X 10"^ < B{B ^ X,7) < 4.5 x 10"^ (95% C.L.). (6.2) 

Employing the leading-order result for B{B Xg'j) (see, e.g., Ref. |]22|), useful bounds can 
be placed on the absolute value of R7, namely 

0.881 < iRj] < 1.321. (6.3) 

B. Integrated observables in B — > K*fi~^ fi~ 

Below we present our predictions (i) for the non-resonant invariant mass spectrum of the 
muon pair, and (ii) for the low-s and high-s region, as defined in Eq. ( |5.4|) . The resonance 
effects are taken into account by employing Eqs. ( p.9[ ) and ( |2.10| ). 

• The branching ratio resulting from the SM operator basis given in Eq. ( p.2|) can be 
conveniently written as 

B = [ao + ai|_RioP + a2|-R7p + Osl-Rgp + a^lmRj + a^lmRg + agRe (-Ryi^g) 

+ arRe R7 + agRe Rg] x 10"^. (6.4) 



^For a review of CP-odd observables, we refer the reader to Ref. [pO[] . 
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TABLE III. Numerical estimate of the coefficients aj = + Pi entering the expressions for the 



integrated branching 


ratio Bi and the averaj 


FB asymmetry (Apg) in B — 


-> K* fj, decay, as 


described in the text. 








Coefficients 






^high 


(ao,/3o) 


(0.03,0.12) 


(0.02,0.08) 


(0.02,0.13) 


(qi,/3i) 


(1.89,8.64) 


(0.62,3.84) 


(0.36,1.76) 


(02, /32) 


(1.05, 1.06) 


(0.95,0.90) 


(0.01,0.04) 


(a3,/?3) 


(1.62,7.34) 


(0.54,3.26) 


(0.30,1.50) 


(a4,P4) 


-(0.07,0.18) 


-(0.01,0.02) 


-(0.02, 0.10) 


(a5,P5) 


(0.23, 0.93) 


(0.01, 0.07) 


(0.11, 0.60) 


{0^6, Pe) 


-(1.52, 3.12) 


-(0.90, 1.59) 


-(0.13,0.52) 


(a7,P7) 


— (0.15, 0.30) 


— (0.12, 0.21) 


(0.02, 0.07) 


(a8,P8) 


(0.33, 1.41) 


(0.18, 0.91) 


-(0.08,0.44) 


Coefficients 








(ao,/?o) 


(0.33,0.33) 


(0.16,0.16) 


-(0.05,0.05) 


(ai,/3i) 


-(1.27,1.42) 


-(0.69,0.77) 


-(0.16,0.18) 


("2,^2) 


(3.29,3.29) 


(0.94,0.94) 


(0.82,0.82) 


(a3,/?3) 


(0.25,0.25) 


(0.01,0.01) 


(0.16,0.16) 



Here we have introduced the coefficients = + (3i which allow us to incorporate 
the effects of new operators [Eq. ( |3.4| )] into our computation of the branching ratio 
by simply replacing 

(a, + ^ (^^f{R^ + K) + l3if{Ri - R[). (6.5) 



where we have defined the quantities R[ analogous to those in Eqs. (|3.1|). Our numerical 
results for the coefficients and (3i are listed in Table |T|. 

Similarly, we parametrize the average FB asymmetry as 

(Afb) = -[Rei?io(ao + alR^ + a2-R9) + a3lmi?io] {~^~^ ' (^•^) 

with the values of Oj tabulated in Table As before, the average FB asymmetry 
in the presence of chirality-flipped operators can be obtained from Eq. ( |6.(j| ) by the 
following replacements: 

(ttj + A)/(-R7, -Rg, -Rio) (^if{.R7 — R'r, R9 — R'g, Rio + R'lo) 

+ PifiRr + R'j, Rg + Ri Rio - R[o). (6.7) 

Our results for the CP asymmetries are as follows: 

(Agp>'°^ = (aoImi?7 + a.lmR,) (^^^ , (6.8) 
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TABLE IV. Values of the coefficients a,; = + /3j for CP-violating observables in i? — > K* ^'^ ^ . 



Coefficients 






(Agp>^°- 


(Agp>^-^ 


(ao,/3o) 


-(0.99,1.76) 


-(0.25,1.03) 


(1.05,1.05) 


(1.62,1.62) 


(ai,/3i) 


(1.20,7.22) 


(1.16,5.98) 







(Agp)''^' = (aoImi?7 + aiIm/?9)(^-gj-j-j, (6.9) 

and 

(Agp>'°" = -aosin0io|i?io|f^Y (6.10) 

<^?P>'''' = -aosin0,o|i?io|(^), (6.11) 

i3 = {B + B)l2 being the CP-averaged branching ratio [see Eqs. (|5.1| ) and ( |5.2| )]. It 
is obvious that (^cp) sensitive to the CP-violating phases 07 and 09, while (^cp) 
also probes the phase 0io. The predictions for the coefficients are reported in Table 

m 



C. Numerical results and predictions 

We first analyse the CP asymmetries in the context of the SM operator basis, Eq. ( ^.21) . 
To determine the implications of new physics for the CP asymmetries, we adopt the following 
procedure. The absolute value of R-j is chosen such that it satisfies the constraints implied 
by the measured b ^ sj rate, Eq. (|6.3| ), while the remaining parameters Rg and Rio are 
required to be consistent with the current experimental upper limit on the non-resonant 
branching ratio B{B^ K*^fi^fi^) given in Eq. ( |6.1|) . To gain predictivity we take {R^l 
to be unity, so that we end up with a set of free parameters comprising (pj, (pg, and \Rg\. 
In fact, imposing the requirement that the non-resonant branching ratio coincides with the 
experimental upper bound, or, alternatively, with the SM prediction, the quantity -Rio is 
computed for any given set of the parameters 07,09, l-Rgj-Q This enables us to consider a 
scenario where the predicted b s'-f fraction coincides with the SM expectation while the 
non-resonant branching ratio oi B ^ K*fi^fi^ may well be accessible in the next round of 
B experiments. 



^As far as new CP-violating phases are concerned, we note that in the context of a specific model 
one has also to take into account the severe constraints on the electric dipole moments of electron 
and neutron. 
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As mentioned above, the value of the strong phase entering the amphtude of the B — >■ 
K*^~^^~ process depends on the Wilson coefficients Ci, . . . , Cg which we have assumed to be 
unaffected by new-physics contributions. Hence, the value of the strong phase sin 6 is fixed, 
and estimated to be —0.07 and —0.51 in the low-s and high-s domain respectively. 

Lastly, using the parametrization for the amplitude given in Eq. ( [5.1| ), and employing the 
integrated expressions above, we determine r and sin0 numerically as a function of 07,09, 
and |i?9|. 

1. CP asymmetry in the partial widths 

Figures ^ and ^ show the parameters r, sin0, and A^p in the low dimuon invariant 
mass region as a function of the phases 07 and 09, taking |i?9| = 1 and requiring that 
= 4.0 X 1O~^.0 It may be noted that the predictions for the average CP asymmetry 
depend very little on the phase 07 and its absolute value is no greater than 1%, regardless of 
the size of the CP- violating phases. If we allow for deviations from |i?9| = 1, the magnitude 
of the CP asymmetry does not change significantly. We may therefore conclude that the 
partial-rate asymmetry, indicative of CP violation, in the low-s region is too small to be 
observable (assuming that the indispensable strong phase does not receive any substantial 
non-standard contributions). 

We now turn our attention to the high-s region. It is clear from the above discussion that 
physics beyond the SM can give rise to sizable CP asymmetries above the ip' resonance where 
we expect to have an appreciable strong phase but also a lower branching ratio. In fact, 
in the high-s region we find the approximate relation r^^^^ ^ 2r^°™, whereas the numerical 
value of the weak phase sin0 is of the same order of magnitude in both the low and the 
high dimuon invariant mass region. This can be seen in Figs. ^ and ^ where we show our 
results for r, sin0, and A^p. For certain values of the phases 07 and 09, the CP asymmetry 
can be of order ±10%. As far as theoretical uncertainties are concerned, we merely remark 
that the numerical estimates for average CP asymmetries are largely independent of the 
parametrizations of form factors, so that the theoretical uncertainty associated with real cc 
intermediate states discussed in Sec. gives by far the largest uncertainty in the predicted 
CP asymmetry. 

Next we consider the case where we allow for higher values of I-R9I. We begin by noting 
that for |_R9| > 1.75 some part of the (07,09) parameter space is already excluded by the 
experimental upper bound on the non- resonant branching ratio, Eq. ( |6.1D . Exploiting the 
fact that A^p depends only weakly on the phase 07, we may take 07 = 4.8 (see Fig. ^. 
Then, setting |i?9| = 1 we find a magnitude of the average CP asymmetry varying from 
—0.05 to 0.12, whereas for |i?9| = 1.75 we predict the range -0.15 ^ (Agp) ^ 0.20. In this 
latter case, numerical values of r between 0.26 and 0.33 are estimated while sin0 can be 
maximal. In particular, for |i?9| = 1.75 the weak phase sin0 is nearly unity if we demand 
the non-resonant branching fraction to be 4.0 x 10~^. 



^The impact of new phases on the partial-rate asymmetry in 5 — > K*l^l decay for the case of 
Ri\ = 1 has also been studied in Ref. p3| . 
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FIG. 4. The parameters r and sm(f> [see Eq. ( ^.31 )] vs and (pg in the low dimuon invariant 
mass region. For simplicity, we have taken |i?7| = \Rg\ = 1 while i?io is chosen to coincide with 
the experimental upper limit on the non-resonant branching ratio, i.e. B^^ = 4.0 x 10~^. 
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FIG. 6. The parameters r and sincj) appearing in ^CP) Eq. ( |5.3D , as a function of cpj and (pg 
for the high dimuon invariant mass region, with |i?7| = \Rg\ = 1. The ratio -Rio is chosen to be 
consistent with the experimental upper hmit on the non-resonant branching ratio given in Eq. ( |6.1[ ). 



0.1- 

high 




27r 271 



FIG. 7. CP- violating average asymmetry (^A^p'^ in the high invariant mass region of the muon 
pair, as a function of (pY and (pg. The parameters \Ri\ (i = 7, 9, 10) have been chosen as in Fig. 0. 
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Let us now discuss a scenario in which we abandon the assumption of having a non- 
resonant branching ratio of 4.0 x 10~^. We focus here on the high-s region as we do not expect 
any significant deviation from our results obtained in the region below the charmonium 
resonances. We first consider the case where the branching ratio is still fixed to the SM 
value of B'^'^ = 1.8 x 10"^. Further, we take 07 = 4.8 and |i?9| = 0.9. (Note that larger 
values of |i?9| are not consistent with the SM branching ratio). As a result, the weak phase 
sin (j) exhibits almost the same 09-dependent behaviour as in the previous case with maximum 
branching ratio presently allowed by experiment. In addition, a smaller value for leads 
to a wider range of r, namely 0.47 ^ r ^ 0.70. As for CP violation, an average asymmetry 
of anything between —0.20 and 0.30 is predicted. 

It is also interesting to analyse the case in which the branching ratio is not fixed to 
any particular value but is still compatible with the experimental results described above. 
Remembering that Rio contributes only to the branching ratio, we set Riq = in order to get 
the highest possible value for the CP asymmetry (see Sec. |V]). Consequently, the interference 
of the terms Rj and Rg now plays an essential role as the branching ratio diminishes for 
certain values of 09, and so r can be unity (which corresponds to the maximal size of the 
CP asymmetry). In this case, the CP asymmetry (^cp) varies considerably and can take 
on any value between —0.5 and 0.4 for B'^^ ranging from 5.0 x 10~^ to 1.6 x 10~^. 



2. CP asymmetry in the angular distribution of /i 

A similar analysis has been carried out for the asymmetry A^p, Eq. ( ^.5|) , which is of 
considerable interest since it probes the phase of Riq. Here again the CP asymmetry in 
the low-s domain is fairly small, typically a few per cent, and we therefore concentrate on 
the high-s region where CP-violating effects are not suppressed by small unitarity phases. 
Results for the ratio Rcp = A^p / sin 0io as a function of 07 and 09 are displayed in Fig. ||. 
Demanding that |i?9| reproduces the experimental upper bound on the non-resonant branch- 
ing ratio leads to the predictions shown in Figure H (a). Moreover, the quantity -Rio obeys the 
constraint \Rio\ ^ 1.9. On the other hand, assuming the SM prediction of B""^^ = 1.8 x 10~^, 
we obtain the results shown in Fig. § (b). In this case, present experimental data on 
B{B^ K*^fi~^fi~) lead to the upper bound \Rio\ ^ 1.2. In both cases, the CP asymmetry 
in the angular distribution of fi~ in B and B decays turns out to be about -10%.^ 

Finally, taking |i?9| = we find that the CP asymmetry can be as large as —25% for 
a rather low |-Rio|, which corresponds to a non-resonant branching fraction of O(10~^). If 
we consider instead a branching ratio of B^^ ~ 10~^, the asymmetry can reach to values of 
-15%. 



^An analysis of such a CP-violating effect in the presence of non-standard Z couplings has recently 
been performed in Ref. |^], which estimates an asymmetry of about 10% in the high dimuon 
invariant mass region. 
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FIG. 8. The quantity i?cp = (^cp) /sin(/)io [cf. Eq. ( |5.1lD ] as a function of ^7 and (j)g 
with lii/l = 1. (a) \Rio\ = 1.9 and |i?9| is chosen such that it coincides with the upper hmit on 
^(5° K*^n+fi-). (b) \Rio\ = 1.2 and the SM branching fraction of 1.8 x 10"^ has been adopted. 



3. A comment on CP violation and additional operators 

As seen in the preceding, within the framework of the SM operator basis it is possible to 
account for the maximum values of the CP-violating asymmetries. Hence our quantitative 
results for CP violation in the decay B K*fi'^fi~ are not affected by the chirality-flipped 
operators [Eq. ( p.4| )] once existing experimental constraints are taken into account. In other 
words, the observation of an appreciable CP asymmetry alone does not provide a test of 
the chirality structure of operators that enter the effective Hamiltonian, and thus is not 
sufficient to disentangle different new-physics scenarios. 

VII. SUMMARY AND CONCLUSIONS 

We have performed a largely model-independent analysis of the exclusive decay B — >■ 
K*l^l~ in the presence of physics transcending the SM. In particular, we have investigated 
the implications of new CP- violating phases for the decay B — ^ K*n~^fi~, and derived 
analytic expressions for the branching ratio, the FB asymmetry, and certain CP-violating 
observables. The formalism presented is applicable to any effective Hamiltonian containing 
the SM operator basis as well as operators with a different chirality structure. We have 
studied in some detail the CP asymmetries in the partial rates and the angular distribution 
of /i~ in 5 — > K*fj,'^fi~ and B K*fi~^fi~ decays, which require the simultaneous presence 
of weak and strong phases. Adopting the SM operator basis and assuming that new-physics 
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contributions to two-body non-leptonic B decays are unlikely to be significant, the CP- 
violating effects in the 2m ^ ^ M^+^- < Mj/^ domain are estimated to be small (up to 
at most a few per cent). Ultimately, this result is a consequence of the smallness of the 
dynamical phase associated with the absorptive part of the penguin diagram; thus, the 
presence of large non-standard CP-violating phases does not necessarily imply sizable CP- 
violating effects in the lower part of the decay spectrum. Even so, studies of CP-violating 
effects in the low-s region will provide a crucial test of the SM, since any indication of CP 
violation would represent new physics. On the other hand, in the high dimuon invariant 
mass region M^/ < M^+^- ^ {Mb — M^*) appreciable CP-violating effects can show up, 
which are consistent with current experimental data on the b ^ sj branching fraction and 
the upper bound on B{B^ K*^jj,'^jj,~). We find that CP asymmetries up to 30% can 
arise for a non-resonant branching ratio of 1.8 x 10~^. That is, new physics gives the same 
rate as in the SM while large CP-violating effects may occur. It should be kept in mind 
that our numerical results for the CP asymmetries in the high dimuon invariant mass region 
are plagued with theoretical uncertainties due mainly to real cc intermediate states, so that 
precise predictions are more difficult in this case. Nevertheless, the CP asymmetries provide a 
particularly useful tool for discovering new physics and their study may gain insight into the 
mechanism of CP violation. Given an asymmetry of 20% and a branching ratio of 3.0 x 10~^ 
in the high-s region, a measurement at 3(j level will necessitate at least 7.5 x 10*^ hb pairs. 



which seems to be achievable in the i?-factory era (see, e.g., Ref. ||25[). In this connection 
it is worth pointing out that the asymmetry A^p, which is a CP- violating effect related to 
the angular distribution of /x~ in B and B decays, has the piquant feature that it does not 
require flavour identification and can be obtained from a measurement of the sum of B and 
B events. 

As far as new operators are concerned, we have argued that in the case of massless 
leptons, i.e. I = e or fi, the dominant contributions may come from operators with a non-SM 
chirality structure. However, since the SM operator basis can accommodate maximum CP 
asymmetries, the inclusion of new operator structures does not affect the main conclusions 
of our analysis, but it is worth considering once sufficient data are accumulated. We are 
thus eagerly awaiting the upcoming B experiments which will provide useful information on 
the short- distance coefficients governing FCNC processes like b —>■ s'-f and b sl^l~ . 

We conclude that the CP asymmetries in the exclusive decay B K*l^l^ can serve as 
an important test of the SM mechanism of CP violation and hence provide a testing ground 
for new physics. A measurement of these asymmetries in forthcoming B experiments would 
signal the presence of non-standard physics and rule out the SM as a primary source of CP 
violation. 
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